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O , Abstract. Semiclassical transport simulation of nucleus- nucleus collisions for the range of 

2 ■ incident energy from about the Fermi energy up to a few hundred MeV per nucleon evidences 

that the maximal excitation energy put into a nuclear system during the early compact stage 
of heavy-ion reaction is a constant fraction of the center-of-mass available energy of the system. 
Analysis of experimental data without presuming reaction mechanism dominating the collision 
^ ' process on the best corroborate the found constancy of energy partition in central heavy-ion 

. reactions. 

o 

O , 1. Introduction 

■ Transformation of the entrance channel longitudinal motion into internal degrees of freedom in 

the course of heavy- ion reaction (HIR) is still awaiting for a satisfactory answer. With increasing 
incident energy Ei^, in particular for central collisions, the reaction mechanism evolves from a 
^ ■ slow, essentially mean-field-transformation fusion and fusion-like processes to a much faster and 

considerably more violent reaction mechanism dominated by elementary nucleon-nucleon (NN) 
collisions. In fusion the entire available energy of the reaction is deposited via thermal excitation, 
whereas at higher energy a considerable fraction of the available energy is deposited into system 
via compression. By increasing E'in one expects that only a fraction of the available energy is 
effectively deposited into the reaction system and becomes dissipated during the reaction course. 
It is commonly admitted that this fraction should monotonically decrease with the increase of 

At Sin from about the Fermi energy Ep to about WOA MeV, the energy transformation is 
determined by those processes which govern heating and compression of a reacting system. In 
this energy range the time scales involved are rather short and are of the order of time which 
reaction partners need to bypass each other [1,2]. Two entrance channel factors play a central 
role in the determination of the dominant reaction mechanism: projectile energy per nucleon 
i^in and reaction geometry (impact parameter and system mass asymmetry). Consequently, the 
course of a HIR is "decided" in the very first instances of a collision [3,4]. In central, the most 
violent collisions the largest fraction of the entrance channel energy is converted into internal 
degrees of freedom. Thus, the central collisions are of our greatest interest. 

We have shown theoretically that an intermediate energy HIR follows a two-stage scenario, a 
prompt first compact-stage and a second after-breakup one [4]. The emission pattern of central 



collisions is characterized by a copious and prompt dynamical emission occurring during the 
compact and prior-to-scission reaction phase [4-6]. This is the main system-cooling component 
and the amount of deposited energy into the compact system linearly increases with the projectile 
energy [7] . These results corroborate conclusion that global characteristics of HIR exit channel 
are determined in the first prompt reaction stage underlying the interest in studying the first 
instances of nuclear collisions. 

In this work we theoretically examine how much of the system energy may be temporarily 
stocked into the reaction system in the form of excitation energy as a function of E'in, system 
size Agys and system mass asymmetry. Four mass symmetric and four mass asymmetric central 
reactions were studied at several energies (see Tab. 1 for a review). Comparison with the 
pertinent results deduced from HIR experiments is presented too. 

2. Theoretical approach 

Simulation was carried out within a semiclassical microscopic transport approach of Boltzmann's 
type using the Landau- Vlasov (LV) model [8] . The highly nonlinear LV equation 

^ + {f,H} = I,Mf) (1) 

is solved by the test-particle method, /(r, p; t) is the one-body density distribution function 
describing the spatio-temporal evolution of the system governed by the effective Hamiltonian H 
consisting of the self-consistent nuclear and Coulomb fields. The Dl-Gl momentum-dependent 
interaction due to Gogny (the incompressibility module -froo=228 MeV and the effective mass 
m*/m=0.67) [9] was used to describe the nuclear mean-field potential. { , } stands for the 
Poisson brackets and /coii is the collision integral. The effects of the Pauli-suppressed two-body 
residual NN collisions are treated on average in the Uchling-Uhlenbeck approximation taking 
the isospin- and energy-dependent free-scattering value for the NN cross section o"nn- Such an 
approach is very successful in reproducing a variety of global experimental dynamical observables 
because they are adequately described by the time evolution of the one-body density. Thus, the 
LV model is especially appropriate for describing the early stages of HIR, when the system is 
hot and compressed. 

The observable studied is the thermal component (heat), i.e. one of the two main intrinsic- 
energy deposition components of the early-reaction-stage energy transformation. Heat is stocked 
into the compact system predominantly by NN collisions which occurs in the overlap zone. This 
is corroborated by the fact that the Pauli principle greatly favors NN collisions involving one 



Table 1. Systems and energies studied for central collisions. 



System 


Incident energ 


y {A MeV) 


40Ar+27Al 


25, 


41, 


53, 


65, 


77, 


99 


36Ar+5SNi 


52, 


74, 


95 








40Ar+i07Ag 


20, 


30, 


40, 


45, 


50, 


75, 100 


40Ar+i97Au 


50, 


75, 


100 






36Aj.+36Ar 


32, 


40, 


52, 


74 






s^Ni+s^Ni 


52, 


74, 


90 








129Xc+120Sn 


25, 


32, 


39, 


45, 


50, 


75, 100 


197Au+197Au 


20, 


30, 


40, 


60, 


80, 


100 



nucleoli from the target and one from the projectile. In the most of cases under study the thne is 
too short for the full relaxation of the pressure tensor and establishment of a global equilibrium 
in momentum space. Therefore, it is more correct to name this component the excitation 
energy E^. Detailed definition of the transformation of the (system) available energy -E^v^ii into 
intrinsic and collective degrees of freedom may be found elsewhere [7, 10, 11]. -^avaii defined 
as the center-of-mass system energy per nucleon -E^vaii ~ ^ (17+^^' '^here Ep/Ap = Ei^ and 
Ap{A'j') is the projectile (target) number of nucleons. 



3. Simulation results 

As an example of the time evolution of excitation energy per nucleon the inset of Fig. 1 shows 
Ex/A for the Au+Au reaction at six energies studied. Within a laps of time of merely 40-75 
fm/c after the contact of colliding nuclei occurring at fm/c the excitation energy per nucleon 
Ey_/A reaches a maximum and then its value decreases almost as rapidly as it increased. As 
expected, the maxima are reached earlier and their height increases and width decreases with 
increasing £"111. The regular and nearly symmetric rise and decrease of E^/A with the reaction 
time is a common behavior for all reactions studied. The observed regularity suggests that 
maxima of E^/A are proportional to the total energy deposited during HIR. 



Figure 1. (Color online.) Simulation results of 
the thermal excitation energy per nucleon E^/A 
for central collisions. 

Inset : Time evolution of E^/A for the Au+Au 
reaction at indicated energies. At each time step 
considered are particles that are bound in large 
fragments, in fact the early compact system. 
Main figure : Excitation maxima (-E'x/^)max as a 
function of system available energy E^:^(^ for mass 
asymmetric (open symbols) and mass symmetric 
(filled symbols) systems studied. The thick grey 
AVAILABLE CM. ENERGY (MeV) line is due to the best linear fit to all data points. 
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We are examining the maximal energy that may be dissipated in HIR. Thus, we take the 
maxima of E^/A which we denote by (Sx/^)max- The value of (-Ex /^) max can readily and 
accurately be extracted from the simulation results. Figure 1 depicts how these maxima depends 
on El:^(i for all studied HIR. Abscissa value is shifted for the threshold, the Coulomb barrier 
energy. With this correction the linear fit over all data points crosses abscissa axis closer to 
the origin of the graph. All data points lie very close to the fit line. One is facing a peculiar 
universal linear rise which is independent of Agys and mass asymmetry in the full and a rather 
large span of Ei^ covered in this study. The linear dependence of (-Ex/^)max on Ei^ is, of course, 
present for each individual system studied. Specifying abscissae in -E^vaii rather than in Ei^ 
merely expresses the mass asymmetric systems on the same footing with those which are mass 
symmetric. 

A universal linear dependence of (-E'x/^)max on -E^vaii ^ ^^^^ nearly exact crossing 

of the origin in Fig. 1 has an important and remarkable consequence: Expressing the value of 
maximal excitation in percentage of the system available energy one obtains that the relative 
fraction of (-Ex/^)max in -^avaii almost constant value except for symmetric systems and 

E'in < Ep. This departure from the constancy occurs because when Ei^ decreases below Ep^ 

^ For mass symmetric systems Ef corresponds to -EavSi ~ 8^ MeV. 
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Figure 2. (Color online.) Top : Ratio of the excitation 
energy and the corresponding E^^^i as a function of this 
same available energy -E^vaii simulation results of 

Fig. 1. 

Middle : Cumulated average number of NN collisions per 
nucleon at the time when (-Ex /^) max is reached. 
Bottom : The average value of the ratio of the potential 
to the effectively realized NN collisions per nucleon at 
the time when (-E'x/^)max is reached. The curves are due 
to a parabolic fit and are merely intended to guide the 
eye. 



the value of the maximum (-E'x/^)max decreases faster than -E^vaii itself is decreasing. This is a 
consequence of an ever slower and slower the early compact system energy transformation as Ei^ 
decreases with an ever more broadened maximum (cf. inset in Fig. 1). Therefore, at these lower 
£'in the maximum (-E'x/^)max is no more proportional on the same manner to the total energy 
deposited in HIR as for Ein > Ep: These (E^/A) 

max cannot be compared with an experimental 
Ex/A of fusion reaction, i.e. of adiabatic-like processes. With this restriction in mind, from 
Fig. 2a) one infers that share of E'x/^ in E^:^^y weekly depends on either reaction system or 
incident energy E^^ and amounts 0.39ib0.03 of -E^vaii- other words, during the early energy 
transformation in HIR the maximal excitation energy that may be deposited in the system is a 
constant which amounts about 40 % of the system available energy. As already discussed, this 
energy dissipation is chiefly owing to NN collisions occurring in the reaction overlap zone. At 
the time when (£'x/^)max takes place the average number of NN collisions per nucleon bound in 
the compact system for each individual system increases linearly and slowly with E^^ (cf. in Fig. 
2b)), whereas Pauli suppression of NN collisions approximately quadratically decreases with Ei^^ 
(cf. in Fig. 2c)). In absolute values the above behaviors are within 20% the same for various 
systems in the energy range studied. With the increase of £'in, the energy transferred by an 
individual NN collision and deposited within the overlap zone increases on the average on such 
a way to remain proportional to E^:^^^. The observed behavior of energy deposit in HIR is a 
consequence of these facts. Let us underline that this constancy of the maximum-of-excitation- 
energy share in available energy is evidenced in the fairly broad range of E^^ (quotient of the 
highest and the lowest E^:^^y covered in the simulation is ~ 9) and it is nearly independent of 
system size (studied is the range of ~ 60 < Asys < ~ 400 nucleons) and mass asymmetry (^p : A^ 
is varied between 1:1 and 1:5). 



4. Comparison with experimental results 

An important question is whether the existing central HIR experimental data support our 
simulation results. Most of the energy put into the system during the early reaction phase 
is released by the emission of particles and light and intermediate mass fragments owing to the 
thermal excitation component E^. At energies below 100 A MeV the compression-decompression 
process contributes a little in the total (kinetic) energy dissipation in HIR [12]. At the instant 
at which the maximum (£'x/^)max is reached a negligible emission occurs. At energies of our 
interest it amounts at most 3-5% of the total system mass [7]. Thus, conjunction of the 
(i?x/^)max with the total (kinetic) energy released in HIR seems to be a natural assumption. 
One must keep in mind, however, that a simulation maximum is reached prior to although very 
close (of the order of ~5-10 fm/c) to the time at which the total momentum distribution becomes 
locally spherical, i.e., the instant at which the local equilibrium has been reached in each part 
of the compact subsystem of bound particles [10]. Nevertheless, the system is far from a global 
equilibrium [7] and comparison with experimental E^/A is not straightforward. In addition, 
one must bear in mind that one should limit the comparison to general trend of experimental 
data, i.e. to the possible constancy of (£'x/^)/-£'avaii ^ ^ function of -E^Xi without seeking to 
reproduce the simulation absolute value. The maximal excitation share of (-E'x/^)max in -E^avaii 
of 40 % is reached during the very first reaction phase and if the same value would be extracted 
from experimental data that could not but be a fortuitous result. Indeed, experimental data is 
registered at an infinite time. Hence, it reflects an integral of the full reaction history. Anyway, 
the simulation maxima (-E'x/^)max should be compared with either the maximal value of E^/A 
obtained in an experiment or with the most probable value of E^/A depending on the nature of 
the distribution. 
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Figure 3. (Color online.) Experimentally 
evaluated total excitation energy per nucleon 
or total dissipated energy per nucleon as a 
function of system available energy. Each 
reaction is represented by its own symbol and 
by color are distinguished different analyses 
of the same reaction. The thick full line 
corresponds to the i^avaii displays the 

upper energy limit which may be reached in 
HIR while the thick dash-dotted line depicts 
30 % of this limit. The only data on total E^/A 
above 100^4 MeV are for the Au-|-Au reaction 
at ^in=150A, 250A and 400^4 MeV [26]. They 
are shown in the inset. The axes aspect ratio 
of both the inset and the main diagram is the 
same so that the slope in both is the same. 



Figure 3 displays a collection of experimental data on ii^xM ^^d total energy dissipated 
in central HIR published in periodics during the last two decades [13-26]. Because energy 



dependence is crucial for our comparison from the figure are dropped all single-energy results. 
Each reaction system is depicted by its symbol while the different measurements of the same 
system are distinguished by color (on line). To avoid of entirely spoiling the figure the error 
bars, typically of 5-15 %, are not displayed. To guide the eye, points belonging to the same 
system and the same analysis are connected and they mostly display close-to-linear dependence 
on ii^avair Unlike the simulation result on (£'x/^)max the experimental data points span a large 
domain of the E^jA vs. plane: The extracted excitations per nucleon lie between one 

third and almost the full accessible system energy i^^vair One may speculate that the different 
approaches used in extracting from experiments the pertinent information on the global energy 
deposition in HIR might be at the origin of these much more scattered results. Indeed, in a 
HIR experiment one does not have a direct access to the excitation energies involved. To obtain 
Ey^jA one needs to reconstruct from detected reaction products the total excitation Ey^ of an 
assumed primary emission source but also the source mass A. There is an evident difficulty 
to restore the break-up stage using exclusively asymptotic experimental information which is 
further obscured by an important role played by primary fragments internal excitation causing 
the in-flight emission. To overcome these uncontrolled issues one has to resort to certain more 
or less justiflcd physical assumptions or/and to use theoretical predictions as a guide for data 
analysis. Anyhow, data analyzed on a same footing seems to fall into much narrower zones of 
the El^-^ vs. plane. 

Linear dependence of E^j A on ^'av™!! is not sufficient to obtain a constancy of its fraction in 
available energy. For this constancy the line passing through data points should also pass close to 
the origin of the -E^v™] vs. Ey^jA plane. As an example in Fig. 4a) are shown results for the Xe-|-Sn 
system which have been extensively studied by the INDRA collaboration. Displayed are five 
analyses of apparently the same data set for E^^ between 25^4 and 50^4 MeV [14, 15, 18, 23, 24]. 
Each analysis have used its own approach in selecting data by centrality and its own philosophy 
reagrding the presumed dominant reaction mechanism used to extract the total excitation Ey^ 
and the primary source mass A. Reported E^j A differ substantially among them. The absolute 
value at the same £'in differs up to 80 %. In addition, some of presumed single-source analyses 
display a rising fraction of E^/A in -Eavaii as -E'in increases [18, 23], other falling fraction as E\^ 
increases [24], whereas the most probable dissipated energy [14] and the total energy loss [15] 
displays a weak if any dependence on E\^. One may argue that various selections of central 
events may reflect different physics and, thus, have a different E^/A. One must admit that 
reported differences, including rising, steady and falling behavior seems to be too big. This 
indicates how delicate these analyses are as well as that at least some of them may contain 
incorrect step(s) in the extraction of the reported values. 

Dissipated energy and total energy loss are the analyses inspired by the kinematical arguments 
and do not require presumption on the dominant reaction mechanism. Their drawback is in their 
applicability to the mass-symmetric systems only. Figure 4b) displays results for all systems 
studied by these two approaches in a fairly broad range of E^in- The total energy loss within the 
error bars gives the same constant value for all four systems studied. These results are rather 
weekly depending on Ei^ and may be considered constant. Another example of cases with the 
constant fraction of E^/A in E'avaii shown in Fig. 4c). Displayed are three single-energy studies 
that carefully accounted for the copious midrapidity emission which occurs during the compact 
and prior-to-scission reaction phase discussed earlier [27-29] as well as the only E^/A result 
reported so far above WOA MeV. Within blast model extracted is the total thermal energy for 
the Au-|-Au reaction from 150^ to 400^ MeV [26]. The abscissae labels above the graph frame 
are relative to the Au-|-Au reaction. These Au-|-Au data have recently been revised [30] but a 
strict linearity of the studied ratio as a function of £'in did not change so that the value of our 
fraction should merely be slightly increased. 
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Figure 4. (Color online.) Ratio of the 
excitation energy and the corresponding -E^vaii 
function of this same available energy 
-^avai'i- Symbols used to distinguish different 
systems are the same as in Fig. 3. 
Top : Five different analysis of the Xe+Sn 
reaction for 25 A < E-^^ < 50 A MeV. 
Middle : Ratio values reported in the analyses 
based on the pure kinematical considerations. 
Bottom : Ratio values reported in analyses 
which thoroughly accounted for the pre- 
equilibrium emission component as well as the 
results on the total thermal energy reported 
above 100^4 MeV and for which the abscissae 
labels above the graph frame are relative to. 



5. Conclusions 

In conclusion, a semiclassical transport model study of the early reaction phase of central heavy- 
ion collisions at intermediate energies has been carried out for a variety of system masses, mass 
asymmetries, and energies below WOA MeV. It has been found that the maxima of the excitation 
energy deposited at this early reaction stage into the reaction system represents a constant 
fraction of about 40 % of the total center-of-mass available energy of the system -Eav™!- In heavy- 
ion experiments extracted total dissipated energy per nucleon and total energy loss deduced 
on kinematical arguments display a similar constancy of their share in the system available 
energy. A similar result may be found in total excitation energy extracted from experimental 
observations under condition that the pre-equilibrium emission is properly accounted for. These 
results indicate that the stopping power of nuclear matter is significant even below the threshold 
of nucleon excitation and that it does not change appreciably when expressed in units of the 
center-of-mass available system energy over a wide range of incident energies. 
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